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situation  is  expected  to  be  worse  for  the  humid  air  cases.  Furthermore, 
neither  experiment  nor  model  addresses  heated  air  channels,  which  are  essen¬ 
tial  to  operational  systems.  Thus,  further  work  in  this  area  is  required. 

The  Medea  experiments  constitute  the  first  channel- tracking  studies  with 
wide  pulse  spacing  (200-500  ps).  Our  measurements  of  the  DC  conductivities  at 
the  optimum  channel  tracking  conditions  are  in  good  agreement  with  the  value 
of  2.5  X  10^  sec”^  predicted  by  simple  theory. 

We  present  as  an  appendix  a  reprint  of  a  paper  resulting  from  our  study 
of  benzene  ion  fragmentation  during  photoionization  with  intense  KrF  laser 
pulses.  We  found  that  at  high  laser  intensities  the  total  ion  yield  was 
reduced  below  the  values  predicted  by  a  rate  equation  model  based  on  cross 
sections  and  intermediate  state  lifetimes  determined  in  our  previous  measure¬ 
ments  at  low  laser  Intensities;  furthermore,  extensive  fragmentation  of  the 
ions  was  observed.  These  results  have  implications  for  laser  channel  tracking 
studies. 
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I  INTRODUCTION 


Objectives 

The  objectives  of  this  contract  were  to  develop  and  apply  diagnostics 
techniques  for  the  propagation  of  charged  particle  beams  in  the  atmosphere, 
including  particularly  the  response  of  the  air  channel  to  passage  of  the 
beam.  These  diagnostics  techniques  fall  into  two  categories:  (a)  measurements 
based  on  air  fluorescence  spectroscopy,  and  (b)  afterglow  conductivity 
measurements  based  on  microwave  techniques.  In  addition,  we  performed  a 
variety  of  auxiliary  experiments  to  support  the  CPB  program,  usually  in  our 
area  of  expertise  in  chemical  physics.  During  this  contract  period,  for 
example,  we  studied  the  fragmentation  of  benzene  during  photolonlzatlon  at 
high  KrF  laser  intensities  in  support  of  the  laser-guiding  technique  used  in 
the  Advanced  Test  Accelerator  (ATA)  at  Lawrence  Livermore  National  Laboratory 
(LLNL). 

Scope  of  Activities 

Our  work  during  this  contract  period  involved  almost  exclusively 
afterglow  conductivity  measurements  performed  in  our  own  laboratory  using 
Febetron  706  electron  beam  excitation  [SDE85a] ,  as  well  as  at  Los  Alamos 
National  Laboratory  (LANL)  using  the  PHERMEX  accelerator  [EDS85]  and  at 
McDonnell-Douglas  Research  Laboratory  (MDRL)  using  the  MEDEA  accelerator 
[SDE85b].  Since  these  studies  have  been  described  in  detail  in  the  references 
just  cited,  this  final  technical  report  will  consist  of  only  a  brief  summary 
and  comparison  of  the  results  from  the  three  sets  of  experiments.  We  also 
include  as  an  Appendix  a  reprint  of  the  publication  on  benzene 
photofragmentation/photoionlzation  resulting  from  our  study  mentioned  above 
[RE85] . 
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Approach 


Our  usual  procedure  is  to  develop  a  diagnostic  technique  in  our 
laboratory  and  to  check  it  out  using  our  Febetron  706  as  an  electron  beam 
source.  This  device  produces  a  3-ns-long  pulse  of  600-'keV  electrons  with  peak 
current  of  about  5  kA  in  a  beam  2  cm  in  diameter.  ’’Beam”  is  a  misnomer,  since 
the  electron  motion  is  generally  governed  by  scattering  so  that  the  current 
density  decreases  rapidly  with  distance  from  the  exit  foil.  The  short  pulses 
and  decreasing  current  with  distance  are  in  contrast  to  accelerators  used  for 
propagation  studies;  hence,  the  electron  beam  excitation  rates  in  Febetron 
studies  are  lower  than  appropriate  to  propagating  beams.  Nevertheless,  these 
studies  establish  a  baseline  for  accelerator  tests,  and  allow  a  more  extensive 
investigation  of  various  parameters  of  the  experiment  than  is  usually  possible 
at  other  laboratories. 

After  establishing  a  technique  in  our  laboratory,  we  use  it  to  perform 
measurements  at  other  laboratories  participating  in  the  CPB  program.  Thus, 
during  this  program  we  made  measurements  on  accelerators  at  LANL  and  MDRL  as 
noted  above. 

The  afterglow  conductivity  measurements  cited  above  were  performed 
primarily  by  using  the  microwave  cavity  perturbation  technique,  which  we  have 
described  in  previous  reports  [Ec83,  SDE85a] .  During  this  contract  period,  we 
also  designed  and  constructed  a  35-GHz  microwave  interferometer  which,  because 
of  time  constraints,  was  used  extensively  only  in  the  PHERMEX  tests  [EDS85] . 
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II  EXPERIMENTAL  TECHNIQUES 


The  propagation  of  an  electromagnetic  wave  through  a  plasma  is  governed 
by  a  scalar  dielectric  constant,  k,  which  is  related  to  the  phase  of  the  wave, 
and  a  scalar  conductivity,  a,  which  is  related  to  the  attenuation  of  the 
wave.  These  parameters  can  be  combined  into  a  complex  dielectric  constant, 


Ic  =  tc  -  j  a/  coSq 


(II. 1) 


or,  alternately,  into  a  complex  conductivity, 


0+  -  l)eQ 


(II. 2) 


The  propagation  can  also  be  expressed  in  terms  of  a  complex  refractive 
index. 


,  -  1/2 
=  M-  -  jX  =  < 


(11.3) 


In  terms  of  the  complex  conductivities,  the  refractive  indices  can  be  written 
as 
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In  a  single-pass  propagation  experiment,  the  phase  change  of  the  probe  beam, 
expressed  as  fringes  per  unit  length,  is 


AS  =  (  Pq  -  ii)  /  \  cm' 


(II. 5) 


whereas  the  attenuation  of  the  probe  beam  is 


a  =  Xiii/c. 


cm 


-1 


(II. 6) 


Thus,  if  an  interferometer  experiment  is  performed  in  which  the  fringe  shift 
and  the  attenuation  are  measured  so  as  to  determine  p  and  x>  respectively,  it 
is  then  possible  to  determine  Oj.  and  from  equation  (II.4). 

An  alternate  approach  that  provides  more  sensitivity  is  the  microwave 
cavity  perturbation  technique,  which  bears  the  same  relationship  to  microwave 
interferometry  as  a  Fabry-Perot  interferometer  does  to  a  Mach-Zehnder 
interferometer  in  optics.  That  is,  the  cavity  creates  a  standing-wave  pattern 
that  Interacts  with  the  plasma  on  a  multipass  rather  than  single-pass  basis. 
The  plasma  causes  changes  in  the  resonant  frequency,  and  the  quality 
factor,  Q,  of  the  cavity  as  given  by  the  Slater  perturbation  formula 


At>) 

0) 


■* 

1  V  •  ^ 


e^o)  f  ^  ^ 

-'v  E  •  E  dv 
o 


(II. 7) 


When  the  perturbation  is  small  enough  that  the  electric  fields  retain  their 
->• 

initial  values,  E  =  E  ,  the  perturbation  formula  can  be  written 


All) 

10 


^  I  o(r)  dv 

e  (0  r  -2,  .  j 
0  J  E  (r)  dv 


-  /  E^(r)  dv 

_ _  plasma _ 

£  (0  c  „2,  .  , 

o  j  E  (r)  dv 

cavity 
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Again,  measurements  of  A5  and  Aw,  which  correspond  to  an  attenuation  and  a 
phase  shift,  allow  determination  of  the  real  and  imaginary  conductivities, 
respectively. 

The  above  equations  outline  the  relationships  between  measured  parameters 
(A5  and  Aw  for  the  cavity,  a  and  AS  for  the  interferometer)  and  the  real  and 
imaginary  conductivities.  The  details  of  the  measurements  and  the  data 
analyses  are  given  in  our  interim  technical  reports  referenced  above.  Now, 
the  complex  conductivity  terms  are  not  directly  useful  (although  the  real 
conductivity  is  approximately  equal  to  the  desired  DC-limit  conductivity  in 
many  cases),  but  they  can  be  used  to  derive  several  plasma  parameters  of 
immediate  interest. 

We  note  first  that  a  plasma  conductivity  has  two  components  —  first, 
from  electrons  interacting  with  positive  ions,  and  second,  from  negative  ions 
interacting  with  positive  ions.  Because  the  negative  ions  move  so  much  more 
slowly  than  electrons,  and  because  the  resonant  frequency  (plasma  frequency) 
of  negative/positive  ions  is  so  much  lower  than  for  electrons  and  positive 
ions,  the  contribution  from  negative  ions  is  usually  neglected.  However,  if 
the  ratio  of  number  density  of  negative  ions  to  niomber  density  of  electrons 
exceeds  the  mass  ratio  of  negative  ions  to  electrons  (which  is  of  order  10^), 
the  contribution  from  negative  ions  will  dominate  the  conductivities.  This 
can  happen  in  the  late  afterglow  of  a  plasma  when  the  electron  attachment  rate 
exceeds  the  electron/ion  recombination  rate,  and  almost  all  the  free  electrons 
are  converted  to  negative  ions.  We  believe  this  is  the  case  in  many  of  our 
measurements. 

The  equations  that  follow  will  be  written  for  electrons,  but  they  are 
equally  valid  for  negative  ions  if  the  appropriate  mass,  velocity,  collision 
frequency,  etc.,  are  used.  Now,  the  complex  conductivity  can  be  written  in 
terms  of  the  electron  mass,  m^;  number  density,  n^;  velocity  distribution 
function,  f(v);  and  momentum  transfer  collision  frequency,  v(v) ,  as 


a 


2 

4it  ^e®  1  df(v)  3  , 

- - - - - j - y  ^y 

3  m  V  +  jw  dv 
e 


(II. 8) 


5 


If  f  and  v(v)  is  known  (for  example,  if  it  can  be  assumed  to  be  Maxwellian  at 
an  electron  temperature  Tg) ,  this  equation  can  be  used  to  calculate  the 
conductivities.  We  are  more  concerned  with  inverting  the  equation,  that  is, 
with  determining  n^  and  f  from  the  measured  o's.  Again,  if  it  is  assumed  that 
f  is  Maxwellian,  the  right-hand  side  of  equation  (II. 8)  can  be  calculated  as  a 
function  of  Tg  for  both  the  real  and  imaginary  terras.  The  ratio  of  these 
terms  can  then  be  formed  and  compared  with  the  experimental  ratios  of  oj./a^. 

A  unique  agreement  between  calculated  and  measured  ratios  validates  the 
asstimptions  and  identifies  an  electron  temperature.  This  temperature  (and 
corresponding  f)  can  then  be  used  to  determine  ng  from  equation  (II. 8)  and  the 
measured  Oj.  or  Finally,  the  real  conductivity  in  the  DC  limit,  can 

be  evaluated  from  equation  (II. 8)  by  letting  oj «  v,  which  leads  to 


,  n  e 
4ii  e 


m 


/t, 


df(v)  3  . 


(II.9) 


Oqq  Is  the  parameter  of  direct  Interest  in  analyzing  the  propagation  of 
charged  particle  beams. 

We  used  this  rather  laborious  procedure  to  analyze  the  results  from  our 
Febetron  experiments  as  reported  in  [SDE85a].  We  also  compared  the  plasma 
properties  calculated  that  way  with  values  calculated  using  the  much  simpler 
constant  collision  frequency  approximation.  As  the  name  suggests,  this 
approximation  consists  of  setting  v(v)  =*  constant  in  equation  (II. 8).  This 
leads  to  the  relations 

Of-  Eo'^^  V  ^*^2  •  (II. 10) 

0)  +  V 

where  the  plasma  frequency,  is  given  by 

0^2  =  n^e^/mgSg  .  (II.  11) 
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With  this  approximation,  the  collision  frequency  can  be  determined  simply  from 
the  ratio  of  real  to  imaginary  conductivity: 


Op/ "  v/o) 


(II,12) 


and  the  plasma  frequency  can  then  be  calculated  by  using  this  value  of  v 
together  with  either  the  measured  real  or  imaginary  conductivity  in  equation 
(II. 10). 

It  should  be  emphasized  that  the  experimental  real  and  imaginary 
conductivities  are  valid  independent  of  the  analysis  techniques  just 
described,  and  independent  of  whether  the  conductivity  is  due  to  electrons  or 
negative  ions.  The  accuracy  of  the  calculated  values  of  collision  frequencies 
and  plasma  frequencies  and,  to  a  lesser  extent,  the  DC-llmit  conductivities, 
varies  with  the  data  analysis  technique  and  the  parameter  regime  (e.g.,  low 
versus  high  pressure) .  Note  that  the  exact  procedure  breaks  down  when  the 
collision  frequency  falls  outside  the  range  expected  for  electrons  (e.g.,  when 
negative  ions  make  the  dominant  contribution) ,  but  that  the  constant  collision 
frequency  approximation  still  yields  valid  collision  frequencies  and  plasma 
frequencies. 

As  indicated  above,  we  analyzed  the  results  of  our  Febetron  experiments 
reported  in  [SDE85a]  using  both  the  exact  procedure  and  the  constant  collision 
frequency  approximation.  The  electron  densities  calculated  using  the 
approximation  were  in  good  agreement  with  the  values  from  the  exact 
calculation  at  low  pressures,  and  progressively  deviated  to  30%  smaller  values 
at  100  torr.  Because  of  its  simplicity,  versatility,  and  generally  acceptable 
accuracy,  we  used  the  constant  collision  frequency  approximation  in  subsequent 
data  analyses. 
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Ill  RESULTS 


Electron  Beam  Parameters 

As  noted  in  the  Introduction,  we  performed  afterglow  conductivity 
measurements  on  the  Febetron  electron^beara  source  at  SRI,  the  PHERMEX 
accelerator  at  LANL,  and  the  Medea  accelerator  at  MDRL.  The  parameters  of 
these  devices  are  summarized  in  Table  1 •  The  PHERMEX  accelerator  produced  a 
sequence  of  10  pulses  spaced  20  ns  apart;  our  results  are  for  the  afterglow 
following  that  entire  pulse  train.  The  Medea  accelerator  can  produce  two 
pulses  spaced  200  ps  or  more  apart.  Most  of  our  results  are  for  single-pulse 
operation,  but  we  did  one  measurement  at  50  torr  air  pressure  with  two  pulses 
spaced  250  ps  apart. 

Summary  of  Real  Conductivity  Histories 

Our  measurements  yield  directly  the  real  and  imaginary  conductivities, 
and  since  the  real  conductivities  are  usually  not  too  different  from  the  De¬ 
limit  conductivities  that  are  of  interest  in  propagation  studies,  we  emphasize 
those  results. 

Figures  1  through  6  present  summaries  of  the  (or,  alternately, 
versus  time  histories  at  each  pressure  tested  for  the  experiments  at  the  three 
facilities  (the  latter  nomenclature  emphasizes  that  these  values  are  for  the 
microwave  frequency  of  these  experiments,  rather  than  the  DC  limit).  The 
first  three  figures  are  for  synthetic  air  (79%  N2,  21%  O2),  whereas  the  last 
three  are  for  laboratory  air,  which  contains  varying  amounts  of  water  vapor. 
These  figures  display  the  general  features  of  afterglow  plasmas.  That  is,  the 
conductivities  decay  approximately  exponentially,  with  the  time  constant  for 
decay  increasing  with  decreasing  pressure  down  to  about  10  torr  and  then 
decreasing  again  with  further  decreases  in  pressure;  the  decay  times  are  also 
considerably  shorter  for  laboratory  air  than  for  synthetic  air. 

In  results  from  the  three  facilities  are  in  reasonably  good  agreement,  as 
can  be  seen  if  the  figures  are  overlaid.  Careful  comparison  shows  that  the 
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FIGURE  3  REAL  CONDUCTIVITY  HISTORIES  FOR  SYNTHETIC  AIR  EXCITED 
BY  THE  MEDEA  ELECTRON  BEAM 
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decays  for  p  <  10  torr  are  longer  in  the  Medea  experiments  than  in  the 
Febetron  experiments.  We  have  shown  [SDE85a]  that  these  low  pressure  decays 
are  governed  by  diffusion,  and  the  presence  of  a  quartz  liner  in  the 
experiments  on  the  Febetron  but  not  on  Medea  makes  diffusion  distances  (and 
therefore  times)  shorter  in  the  former  experiments.  The  differences  between 
laboratory  air  and  synthetic  air  (and  also  some  of  the  minor  discrepancies 
between  results  from  the  different  facilities)  are  due  to  water  vapor,  as  we 
will  demonstrate  below. 

Effects  of  Water  Vapor 

Figures  1  through  6  show  that  conductivity  decays  in  humid  air  are  much 
faster  than  in  dry  air.  We  quantified  this  effect  in  a  series  of  experiments 
at  30  torr  using  Febetron  excitation  [SDE85a].  We  began  by  drying  synthetic 
air  as  well  as  possible  (py^ter  “  0*012  torr),  and  then  systematically 
increasing  the  humidity  in  steps  up  to  Py^ter  “3.3  torr.  The  results  are 
shown  in  Figure  7,  where  we  see  a  factor  of  30  decrease  in  the  decay  time  from 
the  driest  to  wettest  case.  Synthetic  air  as  it  came  from  the  cylinder 
contained  approximately  0.04  torr.  This  extreme  sensitivity  to  moisture 
content  causes  considerable  scatter  in  measurements  of  afterglow  conductivity 
decays,  but  more  important,  it  may  cause  day-to-day  or  site-to-site  variations 
in  the  time  for  the  afterglow  conductivity  to  decay  to  optimum  channel¬ 
tracking  conditions  in  CPB  applications  involving  repetitively  pulsed 
accelerators. 

Electron  Density  Histories 

We  pointed  out  above  that  the  conductivities  measured  in  these 
experiments  contain  contributions  from  both  electrons  and  negative  ions.  We 
generally  cannot  distinguish  between  the  two  sources,  and  our  only  evidence 
for  the  Importance  of  the  ionic  contribution  is  when  the  calculated  momentum- 
transfer  collision  frequency,  v,  is  Inconsistent  with  the  values  that  have 
been  measured  for  electrons  colliding  with  the  molecular  constituents  of 
air.  We  do  find  that  to  be  the  situation  for  many  tests,  usually  involving 
laboratory  air  at  higher  pressures,  and  these  cases  have  been  documented  in 
our  interim  technical  reports. 
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FIGURE  7  REAL  CONDUCTIVITY  HISTORIES  FOR  30  TORR  OF  SYNTHETIC  AIR 
WITH  VARYING  AMOUNTS  OF  WATER  VAPOR 
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The  plasma  frequency  determined  experimentally  can  also  be  a  sum  of 
contributions  from  electrons  and  negative  ions.  However,  as  a  simple  means  of 
comparing  various  results,  we  have  treated  the  plasma  frequency  as  due  only  to 
electrons  and  have  calculated  electron  density  histories  using  equation 
(II. 11).  These  values  are  presented  for  each  set  of  experiments  in  the 
relevant  report.  In  many  cases,  the  histories  exhibit  an  exponential  decay, 
and  we  have  calculated  a  decay  constant  in  each  case.  The  decay  rates  are 
presented  in  Figure  8  for  all  the  air  data.  Since  we  expect  the  dominant 
conductivity  decay  mechanism  in  the  late  afterglow  to  be  electron  attachment 
to  oxygen,  we  have  also  plotted  the  rate  for  this  three-body  process  with  dry 
*’air”  as  the  third  body. 

o 

The  decay  rates  for  pressures  above  30  torr  vary  with  the  same  p"^ 
dependence  as  expected  for  the  electron  attachment  process,  but  the  absolute 
magnitudes  are  6  to  100  times  smaller,  depending  on  the  experiment.  As 
expected  from  the  conductivity  histories,  the  rates  for  dry  air  are  smaller 
than  for  laboratory  air.  Also,  the  rates  from  the  Medea  experiments  are 
slower  than  those  from  the  Febetron  tests. 

Below  10  torr,  ambipolar  diffusion  becomes  a  dominant  plasma  loss 
process,  and  the  decay  rates  approach  a  p"^  dependence.  We  commented  above  on 
geometric  differences  between  the  Febetron  and  Medea  experiments  that  account 
for  their  different  diffuslonal  loss  rates. 

Double-Pulse  Effects 

We  have  made  only  one  measurement  of  conductivity  decays  from  widely 
spaced  pulses  to  date, —  that  on  the  Medea  accelerator  at  50- torr  laboratory 
air  [SDE85b] .  The  results  are  shown  in  Figure  9.  A  careful  comparison  of  the 
decays  from  the  two  pulses  shows  that  the  decay  is  slightly  longer  following 
the  second  pulse  than  following  the  first  pulse.  Most  of  the  difference  can 
be  attributed  to  the  fact  that  the  first  pulse  creates  a  density  channel,  and 
the  reduced  collision  rate  in  that  channel  slows  down  the  kinetic  processes. 

In  addition,  the  elevated  temperature  in  the  channel  changes  the  free 
electron/attached  electron  equilibrium  in  favor  of  free  electrons,  which 
raises  the  conductivity  and  slows  the  decay.  These  effects  will  become  very 
important  for  more  intense  beams  and  a  greater  number  of  sequential  pulses. 
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IV  DISCUSSION 


Status  of  Afterglow  Conductivity  Measurements 

The  general  good  agreement  between  the  three  sets  of  experiments  on 
accelerators  with  different  characteristics  indicates  that  the  characteristics 
of  conductivity  decay  as  a  function  of  pressure  are  known  for  cold  channels. 

We  emphasize  that  these  are  cold  channels,  since  the  highest  temperature  we 
can  document  is  approximately  450  K  in  the  case  of  the  Medea  experiments 
[SDE85b].  We  also  emphasize  that  the  conductivity  histories  are  known  only 
empirically;  attempts  to  model  the  conductivity  histories  using  air  chemistry 
codes  are  limited  to  a  preliminary  effort  by  Ali  [A183] . 

Our  results  to  date  have  established  an  extreme  sensitivity  of  the 
conductivity  decay  rate  to  water  vapor  content  of  the  air.  Furthermore,  the 
decay  rates  for  dry  air  do  not  follow  simple  air  chemistry  models,  and  even 
the  complex  air  chemistry  code  of  Ali  did  not  satisfactorily  predict  the  dry 
air  results  [EDS84].  At  this  time,  the  prospect  of  predicting  the  humid  air 
results  seems  quite  remote,  and  we  have  little  indication  what  effect  channel 
heating  will  have  on  the  conductivity  decays.  This  uncertainty  emphasizes  the 
importance  of  continuing  our  measurements  on  higher  intensity  beams  like  the 
Advanced  Test  Accelerator  at  Lawrence  Livermore  National  Laboratory  in  order 
to  establish  an  optimxjm  pulse  interval  for  channel  tracking  in  realistic 
charged  particle  beam  systems. 

Implications  for  Channel  Tracking 

Simple  theory  [HL84]  indicates  that  the  optimum  DC  conductivity  for 
channel  tracking  is  approximately 


4ti  a/c  =  0.1 


(4.1) 
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where  a  is  the  beam  radius  and  c  is  the  speed  of  light.  For  beam  radii  from 
0.2  to  1  cm,  this  gives  *  (0.2  ^  1.2)  x  10^  sec"^. 

The  only  tracking  experiments  to  date  were  performed  on  the  Medea 
accelerator  at  MDRL  [LR85].  Those  measurements  were  made  over  a  range  of 
pressures  from  1  to  50  torr,  and  optimum  tracking  occurred  at  4  and  8  torr 
when  the  interpulse  spacing  was  200  to  500  |is.  Our  conductivity  measurements 

Q 

indicated  that  was  indeed  near  the  expected  value  of  2.5  x  10  sec  for 
those  pressures  and  afterglow  times  [SDE85b]. 

If  the  decay  rates  observed  to  date  for  cold  beams  are  indicative  of  the 

response  to  more  intense  beams,  the  implication  for  repetitively  pulsed 

accelerators  is  rather  severe.  That  is,  from  Figures  1  or  4,  we  see  that  at 

9  -1 

760  torr  the  conductivity  passes  through  the  10  sec  range  in  «  1  ps. 
Pulsing  an  accelerator  that  quickly  will  surely  present  a  technological 
challenge.  It  would  seem  that  the  conductivity  decay  rate  in  a  hot  channel 
will  have  to  be  very  much  reduced  compared  to  a  cold  channel  to  make  channel 
tracking  a  viable  concept. 

Kinetics  Puzzles 

We  pointed  out  in  discussing  Figure  8  that  the  electron-density  decay 
rates  (really,  the  plasma  frequency  decay  rates)  are  much  longer  than  the 
electron  attachment  decay  rates,  even  though  electron  attachment  should  be  the 
dominant  decay  mechanism  in  the  late  afterglow.  The  differences  are  more 
pronounced  for  dry  air  than  for  wet  air,  and  the  differences  increase  with 
increasing  energy  deposition  rate  in  the  channel.  Similar  effects  were 
observed  in  afterglow  conductivity  decays  from  microwave  breakdown  experiments 
in  air  [SML60]. 

If  the  conductivity  is  in  fact  due  to  electrons,  this  discrepancy  must  be 
due  to  a  combination  of  a  reduction  of  the  attachment  rate  in  the  channel 
environment  and  a  countervailing  detachment  process  that  maintains  the  free 
electron  density.  In  the  first  case,  the  results  of  [MCA72]  show  that  with  an 
air  mixture  of  N2  and  O2  as  the  stabilizing  third  body,  the  electron 
attachment  rate  to  O2  increases  with  electron  temperature  to  0.07  eV  and  then 
decreases,  so  that  the  value  at  1  eV  is  about  1/10 th  of  the  room  temperature 
value.  Since  the  electron  temperatures  deduced  in  our  experiments  are  less 
than  0.2  eV  for  pressures  of  10  torr  and  above  [SDE85a] ,  this  effect  should 
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cause  no  more  than  a  factor  of  two  reduction  of  the  attachment  rate.  In  the 
second  case,  there  has  been  speculation  for  a  long  time  [e.g.,  CPB62]  that 
electron  attachment  in  afterglows  could  be  balanced  by  collisional  detachment 
by  vibrationally  excited  molecules  or  by  metastable  electronically  excited 
molecules  [e.g.,  02(^A)].  Since  the  concentrations  of  such  excited  molecules 
increases  with  electron-beam  excitation  rate,  the  collisional  detachment 
mechanism  would  make  decays  slower  when  the  energy  deposition  rate  was  higher, 
in  accord  with  our  experimental  observations. 

An  alternate  possibility  is  that  the  conductivities  we  measure  are  due  to 
negative  ions  rather  than  to  electrons,  which  we  have  inferred  to  be  the  case 
in  some  of  our  experiments.  Then  the  decay  should  be  due  to  an  ion/ion 
recombination  process  rather  than  electron  attachment.  There  are  two  problems 
with  this  explanation.  First,  the  recombination  can  be  either  of  the 
dissociative  type,  which  is  independent  of  neutral  gas  pressure,  or  of  the 
three-body  type,  which  is  first  order  in  gas  pressure.  Neither  type  should 

o 

create  the  apparent  p  dependence  of  the  decay  rate  exhibited  by  the  results 
of  Figure  8.  Second,  a  recombination  process  leads  to  a  charge  density  decay 
rate  that  has  the  form  n“^  «  t.  This  dependence  can  be  difficult  to 
distinguish  from  an  exponential  decay  when  the  variation  is  only  over  a  small 
range,  but  most  of  our  data  spans  at  least  two  orders  of  magnitude  and  in  most 
cases  the  decays  are  much  closer  to  exponential  than  to  recombination  in 
shape . 

It  will  require  an  extensive  air  chemistry  modeling  study  to  resolve 
these  kinetic  issues  and  to  form  a  basis  for  extrapolation  to  other  channel 
conditions.  The  results  presented  in  our  three  interim  reports  form  a  sound 
benchmark  for  evaluation  of  the  modeling  study. 
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V  CONCLUSIONS 


We  have  summarized  and  reviewed  the  results  from  three  sets  of  afterglow 
conductivity  measurements  that  were  performed  (1)  in  our  laboratory  using 
Febetron  706  excitation,  (2)  at  Los  Alamos  National  Laboratory  using  the 
PHERMEX  accelerator,  and  (3)  at  McDonnell-Douglas  Research  Laboratory  using 
the  Medea  accelerator.  The  results  of  these  studies  were  presented  in  detail 
in  our  interim  technical  reports. 

The  results  of  the  three  studies  are  in  essential  agreement  and  form  an 
empirical  basis  for  predicting  the  afterglow  conductivity  histories  for 
electron  beams  that  do  not  heat  the  air  channel  significantly.  We  find  a 
strong  sensitivity  to  the  humidity  of  the  test  air,  which  necessitates 
consideration  of  that  parameter  in  experiments  and  operational  systems.  The 
dry  air  results  are  in  some  disagreement  with  air  chemistry  model 
calculations,  which  prevents  confident  use  of  those  models  in  their  current 
form  for  predictions.  That  situation  is  expected  to  be  worse  for  the  humid 
air  cases.  Furthermore,  neither  experiment  nor  model  addresses  heated  air 
channels,  which  are  essential  to  operational  systems.  Thus,  further  work  in 
this  area  is  required. 

The  Medea  experiments  constitute  the  first  channel- tracking  studies  with 
wide  pulse  spacing  (200-500  ps).  Our  measurements  of  the  DC  conductivities  at 
the  optimum  channel  tracking  conditions  are  in  good  agreement  with  the  value 

Q 

of  2.5  X  10  sec  predicted  by  simple  theory. 

We  present  as  an  appendix  a  reprint  of  a  paper  resulting  from  our  study 
of  benzene  ion  fragmentation  during  photoionization  with  intense  KrF  laser 
pulses.  We  found  that  at  high  laser  intensities  the  total  ion  yield  was 
reduced  below  the  values  predicted  by  a  rate  equation  model  based  on  cross 
sections  and  intermediate  state  lifetimes  determined  in  our  previous 
measurements  at  low  laser  intensities;  furthermore,  extensive  fragmentation  of 
the  ions  was  observed.  These  results  have  implications  for  laser  channel 
tracking  studies. 
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We  have  carefully  remeasured  the  mass  distribution  of  ions  resulting  from  photoionization  of  benzene  at  248  nm.  The 
fragmentation  occurs  at  lower  intensities  than  previously  indicated,  beginning  below  2  MW/cm-.  At  17  MW/cm  .  the  parent 
ion  constitutes  only  12.5%  of  the  total,  while  at  110  MW/cm^  C^^  constitutes  24%  of  the  yield.  The  functiona  dependence  of 
the  absolute  photoionization  yield  as  a  function  of  laser  power  cannot  be  described  by  a  simp  e  four-level  kinetic  m^el 
successfully  uLd  to  describe  the  absolute  MPI  yield  at  low  power.  Instead,  a  highly  excited  neutral  intermediate  is  postulated 
that  can  either  autoionize  or  photodissociate  upon  absorption  of  an  additional  photon. 


1.  Introduction 

It  is  well  known  that  organic  molecules  produce 
fragment  ions  as  well  as  the  parent  ions  when  they  are 
multiphoton  ionized  at  high  laser  intensities  [1-5]. 
However,  there  is  some  question  about  the  quantita¬ 
tive  accuracy  of  previous  studies  carried  out  using  ex- 
cimer  lasers  because  of  the  difficulty  of  accurately 
specifying  the  laser  intensity  in  the  focal  volume  of 
these  multimode  lasers.  Furthermore,  there  are  no 
quantitative  data  on  absolute  ion  yields  upon  resonant 
two-photon  ionization  at  248  nm  at  higher  intensi¬ 
ties.  Such  data  are  urgently  needed  in  order  to  test  or 
expand  models  for  resonant  MPI  in  polyatomic  mole¬ 
cules.  Because  of  the  importance  of  the  results  to 
studies  of  the  guiding  of  electron  beams  by  preionized 
channels  [6],  we  have  repeated  the  measurements  for 
benzene  photoionized  by  a  KrF  laser,  taking  particu¬ 
lar  care  in  the  determination  of  the  focused  laser  in¬ 
tensity.  In  this  paper  we  report  on  fragmentation  pat¬ 
terns  of  benzene  as  a  function  of  laser  power  at  248 
nm  and  on  absolute  yields  for  resonant  MPI  of  ben¬ 
zene  at  high  laser  power  using  the  low-power  results 
of  Bischel  et  al.  [7]  as  calibration  points. 
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2.  Experimental 

We  used  a  Lambda-Physik  model  EMG  101  excimer 
laser  operated  with  standard  stable  resonator  optics  so 
that  it  produced  180  mJ  pulses  of  14  ns  full  width  at 
half  maximum  in  a  beam  0.7  X  2.2  cm^ .  An  aperture 
of  0.35  cm  diameter  isolated  the  uniform  center  of  the 
beam.  This  aperture  was  focused  with  a  5  cm  focal  length 
lens  unto  a  second  “pinhole”  aperture.  The  calculated 
focal  diameter  at  the  pinhole  was  1.13  X  10“^.  cm. 
However,  transmission  measurements  through  pin¬ 
holes  of  different  diameters  showed  that  the  focal 
diameter  was,  in  fact,  only  marginally  smaller  than 
the  0.025  cm  diameter  pinhole  used  in  all  subsequent 
experiments  to  define  the  focal  diameter  of  the  laser 
beam. 

The  power  dependence  of  the  resonant  two-photon 
ionization  of  benzene  was  studied  by  monitoring  the 
ion  signal  as  a  function  of  the  distance  between  aper¬ 
ture  II  and  focusing  lens  LI  (fig.  1)  with  aperture  II 
being  stationary.  For  all  but  the  highest  power  (2.7 
X  10^  W/cm^),  the  focused  laser  beam  was  vignetted 
by  aperture  II  so  that  its  homogeneous  illumination 
was  achieved  throughout  the  study.  We  estimate  that 
the  maximum  power  transmitted  through  the  pinhole 
was  larger  than  the  nominal  value  quoted  by  about 
10-'20%  due  to  the  fact  that  the  focal  diameter  was 
slightly  smaller  than  the  0.025  mm  diameter  pinhole 
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Fig.  1.  Schematic  of  experimental  apparatus  for  the  measure¬ 
ment  of  the  MPI  fragmentation  of  benzene  at  248  nm.  Ap  I, 
11,  apertures;  LI,  L2,  lenses;  W,  quartz  window;  N,  tungsten 
needle  with  40  tim  diameter  capillary;  G,  grounded  grid;  D, 
field  free  drift  section;  SEM,  secondary  electron  multiplier. 


(aperture  II,  fig.  1).  This  systematic  error,  however, 
applies  only  to  the  highest  power  data  taken  in  each 
run,  because  vignetting  of  the  beam  was  achieved  for 
all  other  (lower)  laser  powers. 

The  pinhole  aperture  was  imaged  with  unit  mag¬ 
nification  into  the  center  of  the  ionization  region 
using  a  30  cm  focal  length  lens.  We  measured  the  laser 
energy  with  a  Scientech  energy  meter  after  the  imaging 
lens  but  before  the  cell  window.  The  maximum  in¬ 
tensity  in  the  focal  volume  based  on  the  diameter  of 
the  pinhole,  the  measured  laser  energy,  and  the  mea¬ 
sured  fwhm  of  the  laser  pulse,  was  110  MW/cm^.  We 
reduced  the  intensity  to  values  as  low  as  4.0  MW/cm^ 
by  moving  the  first  lens  away  from  the  pinhole  so  that 
the  laser  beam  passed  through  a  focus  and  was  diverging 
at  the  pinhole.  We  also  made  measurements  at  2.3 
MW/cm^  using  the  apertured  beam  from  the  laser 
with  no  focusing. 

Benzene  was  introduced  into  the  ionization  region 
through  a  tungsten  capillary  nozzle  of  40  /im  diameter, 
which  was  connected  to  a  reservoir  containing  benzene 
at  lower  than  ambient  temperature,  thus  maintaining 
a  constant  pressure.  The  ionization  chamber  was  main¬ 
tained  at  a  pressure  of  about  10“^  Torr,  so  we  assume 


the  benzene  formed  an  effusive  collisionless  free  jet 
that  diverged  gradually  away  from  the  orifice.  The 
image  of  the  pinhole  was  adjusted  to  produce  the 
maximum  ion  signal  for  each  laser  intensity,  which 
presumably  occurs  when  the  focus  is  just  at  the  exit 
of  the  orifice.  Note  that  the  laser  Rayleigh  range  was 
much  greater  than  the  orifice  diameter,  so  that  the 
intensity  should  have  been  uniform  across  the  width 
of  the  jet. 

As  will  be  discussed  later,  our  irradiation  conditions 
do  not  lead  to  either  saturation  or  photobleaching  of 
the  irradiated  volume,  even  at  the  highest  laser  powers 
achieved.  In  view  of  the  fact  that  we  only  ionize 
10%  of  the  sample  at  the  highest  laser  powers,  we 
took  particular  care  to  align  the  focal  region  of 
the  beam  with  the  jet  to  achieve  maximum  overlap. 
This  was  done  by  reflecting  the  laser  beam  before  it 
entered  the  TOP  mass  spectrometer  through  window 
W  (fig.  1)  and  searching  for  the  position  of  the  image 
of  aperture  II.  The  position  of  this  image  was  brought 
to  overlap  with  the  jet  by  adjusting  the  position  of 
aperture  II  and  focal  lens  L2.  Final  “fine  tuning  ’  was 
achieved  by  maximizing  the  total  ion  signal  through 
laser  beam  steering. 

The  ions  produced  were  measured  using  a  time-of- 
flight  mass  spectrometer  [8]  with  a  2.9  kV  accelerating 
voltage  and  a  secondary  electron  multiplier  detector. 
The  voltage  applied  to  the  SEM  was  2.6  kV  in  most 
experiments,  but  was  reduced  to  2.1  kV  for  experi¬ 
ments  at  the  highest  laser  intensities  in  order  to  avoid 
saturation  effects.  In  addition,  it  was  necessary  to  de¬ 
crease  the  stagnation  presssure  behind  the  nozzle  for 
experiments  at  higher  laser  powers  because  of  space 
charge  effects  on  the  ion  flight  time.  We  changed  the 
benzene  pressure  from  0.8  to  80  Torr  by  cooling  the 
benzene  reservoir  from  room  temperature  to  -30®  C. 
We  assume  that  all  ions  were  detected  with  equal  sen¬ 
sitivity,  regardless  of  mass.  When  the  focal  region  was 
moved  vertically  so  that  the  laser  beam  missed  the  jet. 
the  ion  signal  decreased  essentially  to  zero,  indicating 
the  absence  of  background  ionization. 

The  entire  apparatus  is  shown  schematically  in 
fig.  1. 


3.  Results  and  discussion 

The  distribution  of  ion  masses  measured  at  each 
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lauic  1  A 

Percentage  ion  mass  distribution  resulting  from  248  nm  photoionization  of  benzene  at  selected  laser  intensities  (W/cm^) 


Species 

Mass 

2.3(6) 

3.95(6) 

7.24(6) 

1.64(7) 

3.03(7) 

4,09(7) 

1.11(8) 

C* 

12 

3.25 

9.67 

11.91 

23.89 

CH* 

13 

0.29 

0.47 

0.70 

1.32 

CH2 

14 

0.22 

0.23 

0.26 

0.50 

ch; 

15 

0.15 

C2 

24 

0.65 

3.69 

4.38 

10.25 

CjH* 

25 

0.81 

1.84 

1.93 

3.32 

CjHj 

26 

0.20 

3.52 

3.58 

3.32 

3.33 

4.22 

C2H3 

27 

0.20 

3.52 

2.93 

3.14 

3.33 

4.52 

28 

0.33 

0.28 

0.18 

0.31 

Cl 

36 

6.22 

11.44 

10.30 

12.07 

C3H* 

37 

1.06 

9.13 

9.87 

8.83 

7.23 

C3H5 

38 

1.06 

3.11 

2.34 

2.31 

1,61 

C3H; 

39 

0.89 

1.61 

2.95 

6.22 

5.20 

5.64 

5.22 

C4 

48 

2.96 

2.34 

2.36 

1.61 

C4H* 

49 

0.72 

5.93 

4.35 

4.34 

2.82 

C4h5 

50 

1.32 

6.58 

17.50 

13.06 

12.44 

6.84 

C4H3 

51 

4.13 

12.19 

14.23 

10.04 

10,09 

4,89 

C4H; 

52 

2.69 

3.14 

4.69 

3.56 

3.01 

2,69 

1.20 

CsHj 

63 

0.22 

0.40 

0.94 

0.99 

0.43 

0.70 

Csh; 

76 

2.98 

2.95 

1.01 

0.87 

0.95 

0.40 

CeHs 

77 

13.35 

14.79 

4.45 

4.04 

3.47 

2.01 

CeH? 

78 

89.84 

66.75 

41.36 

11.57 

9.53 

9.15 

5.22 

79 

6.35 

5.95 

3.69 

0.89 

0.87 

0.76 

0.60 

laser  intensity  is  listed  in  table  1,  expressed  as  percent 
of  the  total  yield  at  that  intensity.  The  same  results 
are  shown  in  bar  graph  form  in  fig.  2  and  results  for 
selected  ions  are  plotted  versus  laser  intensity  in  fig.  3 . 

The  results  are  qualitatively  similar  to  those  of 
Reilly  and  Kompa  [2].  We  report  a  maximum  of  22 
different  ions,  compared  to  their  18,  but  the  new  spe¬ 
cies  we  detected  were  present  only  in  low  abundances, 
presumably  below  their  sensitwity  limit.  The  major 
and  important  difference  between  our  results  and 
theirs  is  in  the  intensity  levels  required  to  produce 
given  fragmentation  patterns.  A  visual  comparison  of 
their  mass  distributions  with  ours  indicates  similar 
patterns  when  our  intensities  are  about  one-tenth 
their  stated  values.  As  indicated  above,  we  suspect 
that  the  actual  laser  intensities  in  their  experiment 
were  much  lower  than  they  calculated  because  the 
multimode  excimer  laser  did  not  focus  to  as  small  a 
spot  size  as  expected  from  simple  diffraction  theory. 

120 


Our  own  experience,  described  above,  shows  that  the 
discrepancy  between  the  calculated  and  measured 
beam  diameter  can  exceed  a  factor  of  100. 

In  contrast,  our  results  indicate  somewhat  less  frag¬ 
mentation  at  a  given  laser  intensity  than  do  those  of 
Antonov  et  al.  [3].  They  report  fragmentation  of  ben¬ 
zene  beginning  at  0.05  MW/cm^,  and  their  published 
mass  spectrum  at  2  MW/cm^  shows  similar  levels  of 
fragmentation  as  our  results  at  7  MW/cm^.  They  re¬ 
port  that  CgHg  is  the  dominant  ion  at  least  up  to  20 
MW/cm^,  in  reasonable  agreement  with  our  observa¬ 
tion  that  it  is  dominant  up  to  about  1 5  MW/cm^ .  How¬ 
ever,  their  fragmentation  pattern  is  different  in  detail 
from  those  we  observe  at  any  laser  intensity.  In  par¬ 
ticular,  their  pattern  for  mje  =  76,  77,  and  78  (corre- 
spending  to  CgH^ ,  C5H5 ,  and  CgH^),  does  not  match 
our  results  for  the  same  range  of  laser  intensities.  They 
made  their  measurements  by  reducing  their  laser  beam 
size  using  a  telescope,  and  then  aperturing  the  reduced 
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ION  MASS,  m/e  (amu) 


Fig.  2.  Ion  mass  distribution  resulting  from  photoionization 
of  benzene  at  248  nm  as  a  function  of  laser  power  (14  ns 
fwhm  of  laser  pulse). 

beam.  In  order  to  achieve  high  laser  intensities,  they 
used  up  to  50%  of  the  incident  laser  power.  We  sus¬ 
pect  that  the  beam  was  non-uniform  in  intensity,  so 
that  there  were  “hot  spots”  in  the  reduced  beam.  Their 
fragmentation  patterns  would  then  result  from  the  ad¬ 
dition  of  larger  areas  of  low  intensity  (primarily  re¬ 
sponsible  for  the  parent  ion  signals)  and  localized  areas 
of  high  intensity  (which  would  contribute  the  ma¬ 
jority  of  the  fragment  ion  signals). 

As  suggested  previously  [2],  benzene  photoioniza¬ 
tion  and  fragmentation  at  248  nm  appears  to  occur 
by  way  of  formation  of  the  parent  ion  (CgHg),  which 
subsequently  absorbs  more  photons  and  decomposes 
stepwise,  even  as  far  as  C'*’  at  the  highest  intensities. 
The  yield  of  the  parent  ion  decreases  steadily  over  the 
entire  intensity  range.  In  addition,  and  C^H^, 
which  first  appear  at  4  MW/cm^,  decrease  in  constant 


Fig.  3.  Ion  mass  distribution  for  the  most  important  fragment 
groups  resulting  from  photoionization  of  benzene  at  248  nm 
as  a  function  of  laser  intensity. 

proportion  to  at  higher  intensities,  as  shown  in 
fig.  2.  The  similarity  of  appearance  intensities,  together 
with  their  near  equality  of  appearance  energies,  sug¬ 
gests  that  these  two  ions  are  both  formed  indepen¬ 
dently  from  the  parent  ion.  Furthermore,  their  con¬ 
stant  proportionality  to  the  parent  ion  suggests  that 
these  two  product  ions  undergo  photodissociation  at 
the  same  rate  as  the  parent  ion.  As  can  be  seen  in 
figs.  2  and  3,  the  details  of  the  photofragmentation 
spectrum  of  are  a  complex  function  of  the  laser 
intensity.  It  appears  that  elimination  of  CH^j  {n  =  1-3) 
to  yield  C5H3  is  not  a  favored  process  over  the  entire 
intensity  range.  Also,  the  elimination  of  C2H„  (^  “ 
0—4)  seems  to  be  favored  over  the  extrusion  of  C3H„ 
(«  =  0—3)  at  intensities  below  10^  W/cm^.  The  situa¬ 
tion  is  complicated  by  the  fact  that  the  precursors 
(ionic  or  neutral)  for  the  fragment  ions  are  not  present¬ 
ly  known.  However,  the  general  trend  can  be  noted, 
that  the  appearance  energy  for  a  fragment  ion  is  in¬ 
versely  proportional  to  its  carbon  content  (with  the 
exception  of  C5H3). 

We  have  added  the  contributions  of  ions  of  all 
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Fig.  4,  Absolute  ion  yield  (mass  integrated)  in  resonant  two- 
photon  ionization  of  benzene  at  248  nm  as  a  function  of  laser 
power  /  (14  ns  fwhm  of  laser  pulse).  •,  •  ...unfocused  beam, 
80  Torr  stagnation  pressure,  o,  a  ...focused,  0.8  Torr  stagna¬ 
tion  pressure,  □  ..Tocused,  25  Torr  stagnation  pressure. 


masses  at  each  laser  intensity  to  determine  the  de¬ 
pendence  of  total  photoionization  yield  on  laser  power 
or  intensity.  The  results  are  compared  in  fig.  4  with 
the  predicted  variation  based  on  the  molecular  con¬ 
stants  for  benzene  measured  recently  by  Bischel  et  al. 
[7  ] ,  using  our  laser  pulse  shape . 

In  their  work,  absolute  ionization  yields  for  the 
two-photon  resonant  excitation  of  benzene  were  mea¬ 
sured  at  laser  intensities  low  enough  that  no  fragmen¬ 
tation  occurred.  The  ion  cunent  was  found  to  be  strictly 
second  order  in  laser  intensity  up  to  laser  fluence  of 
8  mJ/cm^ ,  and  the  results  were  fitted  to  a  four-level 
rate  equation  model  of  benzene.  The  assumption  was 
made  that  photoionization  only  takes  in  benzene  re¬ 
sulting  in  benzene  molecular  ion,  ,  and  therefore 

ionization  of  neutral  fra^ents  originating  from  the 
photodissociation  of  C5Hg  or  C5Hg  would  not  con¬ 
tribute  to  the  total  observed  ion  current.  This  assump¬ 
tion  also  leads  to  the  staternent  that  all  fragment  ions 
must  have  their  origin  in  ,  an  assumption  that 
may  not  hoM  in  view  of  the  low  ionization  potential  of 
many  open-shell  hydrocarbon  species.  In  fact,  Bischel 
et  al.  were  able  to  explain  their  data  on  the  basis  of 
the  above  assumption  (no  fragmentation  under  their 
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experimental  condition) 
model: 

^A|g  +  /zJ^ — ^  ^®2u  ’ 

,  ^21  + 
1B2u— 


and  the  following  four-level 

(1) 

(2) 

(3) 


where  is  the  electronic  ground  state,  ^B2u  is  the 
resonant  intermediate  state,  ^Bj^  is  the  non-photo- 
ionizing  lowest  triplet  state,  and  is  the  ionic 
ground  state  of  benzene;  Ojo  and  031  are  the  corre¬ 
sponding  absorption  cross  sections,  and  Q  is  the  effec¬ 
tive  rate  of  deactivation  of  the  resonant  intermediate 
state  with  radiative  (fluorescence)  and  non-radiative 
(intersystem  crossing)  components.  With  the  aid  of 
known  values  for  of  3.7  X  10“^^  cm^  and  Q  of 
2.2  X  10^  an  ionization  cross  section  of  4.4X10“^^ 
cm^  was  found  for  031 » together  with  a  small  non- 
resonant  two-photon  cross  section  of  <  0.5  X  10^^^ 
cm^/W  [7]. 

The  solid  line  in  fig.  4  describes  the  absolute  ioni¬ 
zation  yields  using  the  above  molecular  parameters  for 
benzene  and  our  laser  beam  temporal  profile.  The  low- 
intensity  yields  (black  symbols  in  fig.  4)  were  then 
fitted  to  the  calculated  curve,  because  the  absolute 
scale  for  the  experimental  points  is  arbitrary.  The 
low-intensity  data  fit  a  second-order  dependence  on 
laser  intensity  rather  well  up  to  the  point  where  frag¬ 
mentation  begins  to  become  important  (5%  of  total 
ion  current  at  2.5  X  10^  W/cm^  is  distributed  to  lower- 
weight-fragment  ions,  see  table  1).  The  high-intensity 
mass  integrated  ion  yields  (open  symbols  in  fig.  4) 
were  fitted  such  that  the  lowest  two  intensity  data 
points  coincided  with  the  prediction  from  our  rate 
equation  model.  The  experimental  results  at  intensi¬ 
ties  exceeding  7.5  X  10^  W/cm^  fall  well  below  the 
prediction,  which  suggests  either  that  an  additional 
loss  process  “turns  on”  at  higher  laser  power,  or  that 
space-charge  effects  reduce  the  collection  efficiency 
for  all  ions  at  higher  ion  densities.  However,  the  rela¬ 
tive  ionization  yield  for  a  variation  of  the  benzene 
backing  pressure  of  as  much  as  a  factor  of  100  was 
found  to  be  unchanged,  which  we  take  as  an  indica¬ 
tion  that  space  charge  effects  or  recombination  ki¬ 
netics  were  not  responsible  for  the  depression  of  the 
MPI  yield  at  higher  laser  power. 
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At  this  point  our  results  lead  us  to  postulate  a  neu¬ 
tral  precursor  state  that  can  either  undergo  autoioniza¬ 
tion  or  absorb  an  additional  photon  to  decay  into 
neutral  fragments.  The  depression  of  the  MPI  yield  at 
higher  laser  intensities  cannot  be  due  to  saturation  ef¬ 
fects  because  they  have  already  been  taken  into  account 
in  the  rate  equation  model.  Furthermore,  saturation 
effects  are  unimportant  in  the  present  case  because 
Oji  >  lOojQ.  This  state  is  thought  of  as  a  superex- 
cited  state  of  benzene  at  an  energy  of  twice  the  photon 
energy  above  ground  state  and  is  most  likely  a  mo¬ 
lecular  Rydberg  state  of  benzene.  Using  this  neutral 
intermediate  the  resonant  MPI  scheme  is  given  by 


M  +  M*  ,  (4) 

+  (5) 

(6) 

+  (7) 

M**-^M+  +  e,  (8) 


where  M  represents  C5H5 ,  M*  the  one-photon  reso¬ 
nant  intermediate,  and  M**  the  superexcited  state  of 
benzene. 


4.  Conclusions 

We  have  reported  careful  quantitative  measurements 
of  the  distribution  of  ions  resulting  from  photoioniza¬ 
tion  of  benzene  with  a  KrF  laser  at  248  nm.  The  frag¬ 
mentation  is  extensive,  as  reported  previously  by  Reilly 
and  Kompa,  but  our  measurements  indicate  that  it 
occurs  at  about  one-tenth  the  laser  intensity  that  they 
reported.  Our  results  are  in  reasonable  agreement  with 
those  of  Antonov,  Letokhov  and  Shibanov,  although 


their  intensities  for  comparable  fragmentation  were 
even  lower  than  ours.  We  have  also  reported  absolute 
MPI  yields  at  high  laser  powers  which  led  us  to  postu¬ 
late  a  neutral  superexcited  molecular  state  of  benzene. 
The  presence  of  this  superexcited  state  is  thought  to 
be  responsible  for  the  depression  of  the  absolute  MPI 
yield  with  respect  to  the  predicted  yield  for  laser 
powers  exceeding  7.5  X  10^  W/cm^. 
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